Scaling in the Lattice Gas Model 
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A good quality scaling of the cluster size distributions is obtained for the Lattice Gas Model using 
the Fisher's ansatz for the scaling function. This scaling identifies a pseudo-critical line in the phase 
diagram of the model that spans the whole (subcritical to supercritical) density range. The indepen- 
dent cluster hypothesis of the Fisher approach is shown to describe correctly the thermodynamics 
of the lattice only far away from the critical point. 



PACS numbers: 24.10.Pa,64.60.Fr,68.35.Rh 
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Since the first heavy ion experiments multifragmen- 
tation has been tentatively connected to a critical phe- 
nomenon p]. The recent determination of a consistent 
set of critical exponents in different multifragmentation 
data 0, [H tends to confirm this hypothesis even if the 
finite size corrections to scaling are largely unknown. On 
the other side the experimental observation of a flatten- 
ing of the caloric curve [0, ^ and the measurement of 
a negative heat capacity points towards a first order 
phase transition as it is also suggested by the thermody- 
namics of statistical multifragmentation models j?J . The 
debate on the order of the transition has been further an- 
imated by a very recent analysis of the Isis data || which 
shows a high quality scaling of the fragment size distribu- 
tion over a wide range of charges and deposited energies 
with an ansatz for the scaling function taken from the 
Fisher droplet model B which approximates a real fluid 
as an ideal gas of clusters. The critical temperature ex- 
tracted from the Fisher analysis is identified as the tem- 
perature of the thermodynamical critical point and the 
whole coexistence line of the nuclear phase diagram is re- 
constructed under the hypothesis that the Fisher model 
gives a good description of the multifragmentation phe- 
nomenon ||]. 

In this paper we apply the Fisher scaling method of 
ref. H to the Lattice Gas model, which is a well known 
paradigm of first as well as second order phase transitions 
JlO| |, the canonical version of this model being isomor- 
phous to the Ising model at fixed magnetization. We will 
show that the observation of Fisher scaling does not allow 
to determine the location of the critical point. Critical 
behaviors are observed on a line at supercritical density 
PI which extends inside the coexistence region [D| 
due to finite size effects. We suggest that the inadequacy 
of the Fisher model to describe the lattice gas phase dia- 
gram comes from the Fisher independent cluster hypoth- 
esis which should be correct only at low densities and high 
temperatures ; indeed we demonstrate that the thermo- 
dynamics reconstructed from the Fisher partition sum 
coincides with the exact one only in the low density and 
high temperature regime. 

In our implementation of the lattice gas model jl0| the 
N sites of a cubic lattice are characterized by an occupa- 
tion number n, which is defined as = 0(1) for a vacancy 
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FIG. 1: RG scaling from eq§(figs. la,lc) and Fisher scal- 
ing from eq.^ (figs. lb, Id) of the cluster size distribution in 
a 8X8X8 cubic lattice at the critical density for temperatures 
0.36 < T/c < 3.6 and cluster sizes 4 < A < 30. The critical 
parameters corresponding to each figure are indicated. 



(nucleon). Particles occupying nearest neighboring sites 
interact with a constant coupling c. The relative particle 
density p/po is defined as the number of occupied sites 
divided by the total number of sites and is linked to the 
magnetization of the Ising model by p/ po = 2m — 1. Ob- 
servables expectation values are evaluated in the canoni- 
cal ensemble sampled through a standard Metropolis al- 
gorithm Q. The use of the canonical constraint allows 
a direct exploration of the coexistence region || . The 
backbending of the chemical potential as a function of 
the particle number and of the pressure as a function of 
the density allows an unambiguous definition of the phase 
diagram even for a finite system 1 12 . Fragments are de- 
fined within the standard Coniglio-Klein bond breaking 
probability between occupied neighbor sites. For all tech- 
nical details, see refs.|l^, 

Renormalization group arguments lead to the expecta- 
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tion that in the vicinity of the critical point (T c , p c ) the 
size distribution scales as Jl5| 

n{A,T) = qA- T f {eA°) (1) 

where e measures the distance from the critical point 
e = (T c — T)/T c , q is a normalization constant, / is the 
scaling function and r, a are critical exponents. We will 
refer to eq.(||) in the following as to RG scaling. The RG 
scaling analysis performed on a 8X8X8 lattice with peri- 
odic boundary conditions at the critical density p c = 1/2 
is displayed in figure lajl2|, [l4|. A good scaling behav- 
ior is observed for all temperatures 0.36 < T/c < 3.6 
and all cluster sizes 4 < A < 30. The critical exponents 
t = 2.2,cr = 0.61 are close to the expected values of 
the liquid-gas universality class r = 2.2, a = 0.64 and 
the critical temperature T c = 1.16c is in good agreement 
both with the temperature at the thermodynamical crit- 
ical point jlj] T^ h — 1.22c and with the expected critical 
temperature at the thermodynamical limit = 1.12c. 
Indeed finite size corrections to scaling have been eval- 
uated HH and found to be small. The method used to 
extract the critical parameters is discussed in detail in 
refs. 0,0. 

In the Fisher droplet model |J the vapor coexisting 
with a liquid in the mixed phase of a liquid-gas phase 
transition is schematized as an ideal gas of clusters. A 
similar scaling around the critical point is supposed by 
this model but a different form is suggested for the scaling 
function 

n(A,T) = q A-Texp( AA »- T CoeA(T ) (2) 

Here A/j, represents the difference in chemical potential 
between the two phases, and cq is the surface energy coef- 
ficient. Since both A/z and cq can be in principle temper- 
ature dependent, we have parametrized these quantities 
as poly nomials of order 4 and 1 respectively following 
rcf. |16| ; the normalization q has been taken as in the 
infinite system JlTf . The critical parameters obtained 
from the best )C".t as well as the scaled distributions 
are shown in figure lb. Only temperatures lower than 
the maximum production temperature for each size have 
been used in the fit (black dots in figure 1). The scal- 
ing is violated only for higher temperatures (grey dots) 
consistently with the Fisher approach which modelizes 
only the vapor coexisting with a liquid, i.e. is relevant 
for temperatures T < T c . It is surprising that two such 
different ansatz for the scaling function lead to a compa- 
rable quality for the scaling of the size distributions and 
to coherent and close values for the critical exponents; 
this remarkable result confirms the wide universal valid- 
ity of generic thermal scalings |l8| . The main ambiguity 
concerns the critical temperature which comes out about 
20% higher with the Fisher technique. This difference is 
not a broadening effect due to finite sizes, a well defined 
critical point being replaced in small systems by a wide 
spread critical region. As a matter of fact, if the RG 
critical parameters of figure la are implemented in the 
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FIG. 2: Fig. 2a: Fisher scaling as in fig. lb but at different den- 
sities. Fig. 2b: thermodynamical coexistence line (full line) 
and region of critical partitions (dashed lines) from ref.[E3]]. 
Grey lines: coexistence line reconstructed from fragment par- 
titions via eq.(3). 



Fisher analysis and a reduced 7 parameters fit is done 
with the ansatz (0) the scaling is clearly violated even at 
low temperatures (figure Id) and the same thing is true if 
the Fisher scaling parameters of figure lb are inserted in 
the RG ansatz eq.(Q) as shown by figure lc. This means 
that the two scaling ansatz are not equivalent and the 
good quality of the scaling is insufficient to prove the 
adaptation of the model to the data. Therefore before 
giving a physical meaning to the precise value of the ex- 
tracted critical temperature one should a priori know if 
the chosen scaling ansatz is consistent with the system 
under study. On the other side the critical exponents 
seem to be very robust and depend only very slightly on 
the scaling hypothesis. 

The analysis shown in figure 1 was performed at a con- 
stant density equal to the critical density.In the case of 
nuclear collisions it is not obvious that multifragmen- 
tation occurs systematically at the critical density. In 
particular the good adequacy between statistical models 
and data favours a lower value for the freeze out den- 
sity JtJ and similar information come from interferometry 
measurements |lj|. On the other side recent calculations 
in the framework of classical molecular dynamics [Tl} | pro- 
pose an early fragmentation at supercritical density. To 
understand the effect of volume we have performed dif- 
ferent lattice gas calculations at different densities. As 
shown in the upper part of figure 2, a very good scal- 
ing is observed for all subcritical as well as supercritical 
densities. In all cases the values of the critical exponents 
are comparable but the critical temperature is a mono- 
tonically increasing function of the density. To visualize 
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all the results on the same picture a constant horizon- 
tal shift C(p) is given to each scaled distribution. The 
critical temperatures obtained for each density are repre- 
sented by the black symbols in the lower part of figure 2. 
In this figure the full line gives the coexistence line of the 
model calculated in a precise way from the derivatives 
of the canonical partition sum [Q. The locus of criti- 
cality lies approximately over a line which passes close 
to the thermodynamical critical point (open dot) but ex- 
tends further at supercritical (Kertesz lineal)) as well 
as subcritical densities inside the coexistence region. A 
qualitatively similar behavior has been already observed 
with the RG analysis (dashed line in figure 2b) |L2| and 
has been interpreted in terms of finite size effects. This 
result implies that in the framework of the lattice gas 
model the observation of Fisher scaling and more gener- 
ally of a critical behavior does not allow to localize the 
critical point and is compatible also with fragment for- 
mation at low density inside the coexistence region. 

A first order phase transition in a finite system cor- 
responds to a concavity anomaly in the free energy 
F = TlnZ which in turn leads to a backbending of the 
canonical chemical potential fi = f + pd p f where / is the 
free energy per particle. The coexistence line in figure 
2b corresponds to the equality of the chemical poten- 
tials A/i = on the liquid and gas branch defined by a 
Maxwell construction jl2|, ^0|. If the Fisher model is a 
good approximation to the Lattice Gas physics it should 
be therefore possible to reconstruct the vapor side of the 
coexistence line pcl (T) directly from the fragment yields 
as proposed in ref.|| as 

Pcl = n ^ T)Aexp(-AAp) (3) 

A 

where the sum extends over all fragments but the biggest. 
The resulting curves are given by the grey lines in fig- 
ure 2b for the four different densities shown in figure 
2a. The end point of the lines giving by construction 
the total density of the system and the critical tempera- 
ture extracted by the Fisher fit, these lines are obviously 
meaningless if the system is fragmenting at a density dif- 
ferent from the critical density p c = 1/2. However even 
at p = p c when the thermodynamical critical point is 
included in the data set the reconstruction of the coex- 
istence line is very poor. In particular the curvature of 
this line at the reconstructed critical temperature corre- 
sponds to an exponent f3 — 0.84 which strongly deviates 
from the expected (3 = 0.31 exponent given by the ther- 
modynamical coexistence line (full line in fig. 2b) which 
would be consistent with the liquid gas universality class 
(/3 = 0.33)and with the critical exponents extracted from 
the clusters. 

This means that despite the magnificent scaling shown 
by figure 2a the physics of the Fisher droplet model does 
not correspond to the Lattice gas. 

This may look surprising since the main hypotheses 
of the Fisher model are shared by the Lattice : clus- 
ters are essentially defined by a volume and a surface 
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FIG. 3: normalized pressure versus temperature at two differ- 
ent densities from the exact canonical partition sum (full dots) 
and from the ideal gas approximation eq.(4) (open dots). 



contribution; they exhibit a critical behavior at the ther- 
modynamical critical point; the statistical weight of a 
given configuration is given by a Boltzmann factor. The 
spectacular collapse of all the cluster distributions on the 
single curve of figure 2 indeed indicates that the Fisher 
ansatz gives a good prediction of an important part of the 
physics of the Lattice, i.e. the inclusive yields. However 
it may be interesting to remark that a somewhat differ- 
ent value for j3 can be obtained by changing the order 
of the polynomial assumed for A/i without any sizeable 
change in the quality of the scaling. This suggests that 
the informations contained in the inclusive yields may 
be insufficient to pin down the thermodynamics of the 
system. 

The two models strongly differ on one point: if in the 
Fisher picture fragments constitute an ideal vapor of non 
interacting composite particles (the individual produc- 
tion probabilities are factorized) in all microscopic mod- 
els as the Lattice Gas, interactions among fragments are 
naturally taken into account through the volume they 
exclude and through the surface coupling between neigh- 
boring fragments. One may therefore wonder if these in- 
teractions, which seem to affect in a non crucial way the 
inclusive yields, may induce important differences in the 
thermodynamics. The importance of these effects can be 
studied by testing the deviation of the Lattice equation of 
state from the ideal gas hypothesis of the Fisher model. 
If fragments can be modelized as an ideal classical gas, 
in a constant volume transformation the pressure can be 
deduced directly from the fragments yields [|| 

P_ = EW^T) T 

Pc £ A n(A,T c )T c y> 
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where the sum extends over the vapor phase (all frag- 
ments but the biggest) and T c — T c (p) is the temperature 
obtained from the Fisher fit for each density (black dots 
in figure 2b). This pressure is shown by the open dots 
in figure 3 at two different densities and can be com- 
pared to the exact pressure of the model from ref.flifl 
p = p 2 d p f (black dots). Not surprisingly, the gas of clus- 
ters behaves as an ideal gas only at low density and high 
temperature. The attractive interaction among fragment 
surfaces causes the pressure to become negative at low 
temperatures (the system is bound) while an extra pres- 
sure comes from the excluded volume interaction at high 
density. 

In conclusion in this paper we have analyzed the frag- 
ment size distributions issued of the canonical implemen- 
tation of the Lattice Gas model by means of Fisher scal- 
ing. A very good scaling is observed at subcritical as 
well as supercritical densities with values for the critical 
exponents compatible (within finite size effects) with the 
universality class of the model. This implies that the ob- 



servation of scaling does not allow to infer the position of 
the critical point and is compatible with a fragmentation 
inside the coexistence region of a first order phase tran- 
sition. Knowing that the scaling function of the model 
is sensibly different from the Fisher exponential ansatz, 
Fisher scaling appears as a very powerful tool to extract 
critical exponents in a way which is essentially indepen- 
dent of the assumptions made on the detailed shape of 
the scaling function. On the other side the reliability 
of thermodynamical quantities extracted from the Fisher 
analysis (coexistence curve, critical temperature, satu- 
rated pressure..) for any set of experimental or simulated 
data depends critically on the possibility of approximat- 
ing the fragment partitions as an ideal noninteracting 
gas. In the case of the Lattice Gas model this hypothesis 
is verified only for densities much lower and temperatures 
much higher than the ones of the critical point. 
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